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ABSTRACT 



The new cosmic microwave background (CMB) temperature maps from Planck provide the highest-quality full- sky view of the surface of last 
scattering available to date. This allows us to detect possible departures from the standard model of a globally homogeneous and isotropic cosmol- 
ogy on the largest scales. We search for correlations induced by a possible non-trivial topology with a fundamental domain intersecting, or nearly 
CO intersecting, the last scattering surface (at comoving distance ;^rec), both via a direct search for matched circular patterns at the intersections and 
by an optimal likelihood search for specific topologies. We consider flat spaces with cubic toroidal (T3), equal-sided chimney (T2) and slab (Tl) 
CO topologies, three multi-connected spaces of constant positive curvature (dodecahedral, truncated cube and octahedral) and two compact negative- 
^-H curvature spaces. These searches yield no detection of the compact topology with the scale below the diameter of the last scattering surface. For 
!• most compact topologies studied the likelihood maximized over the orientation of the space relative to the observed map shows some preference 
• multi-connected models just larger than the diameter of the last scattering surface. Since this eff'ect is also present in simulated realizations of 

k>< isotropic maps, we interpret it as the inevitable alignment of mild anisotropic correlations with chance features in a single sky realization; such a 
feature can also be present, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius % 
^ of the largest sphere inscribed in topological domain (at log-likelihood-ratio AlnX > -5 relative to a simply-connected flat Planck best-fit model) 
are: in a flat Universe, > Q.92xxqc for the T3 cubic torus; Hi > ^.Hxxqc for the T2 chimney; Hi > 0.50;^rec for the Tl slab; and in a positively 
curved Universe, Hi > 1.03;^rec for the dodecahedral space; Hi > l.^Xrcc for the truncated cube; and Hi > 0.89;^rec for the octahedral space. The 
limit for the T3 cubic torus from the matched-circles search is, consistently, Hi > 0.9Axxqc at 99 % confidence level. 

We also perform a Bayesian search for an anisotropic global Bianchi VII/^ geometry. In the non-physical setting where the Bianchi cosmology 
is decoupled from the standard cosmology, Planck data favour the inclusion of a Bianchi component with a Bayes factor of at least 1.5 units of 
log-evidence. Indeed, the Bianchi pattern is quite efficient at accounting for some of the large-scale anomalies found in Planck data. However, the 
cosmological parameters that generate this pattern are in strong disagreement with those found from CMB anisotropy data alone. In the physically 
motivated setting where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no 
evidence for a Bianchi VII/j cosmology and constrain the vorticity of such models to (a)/H)o < 8.1 x 10~^^ (95 % confidence level). 

Key words, cosmology: observations — cosmic background radiation — cosmological parameters — Gravitation — Methods: data analysis — 
Methods: statistical 
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1. Introduction 

This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration I 2013), de- 
scribes the use of Planck data to limit departures from the global 
isotropy and homogeneity of spacetime. We will use Planck'^ 
measurements of the cosmic microwave background (CMB) to 
assess the properties of anisotropic geometries (i.e., Bianchi 
models) and non-trivial topologies (e.g., the torus). The sim- 
plest models of spacetime are globally isotropic and simply con- 
nected. Although both are supported by both local observations 
and previous CMB observations, without a fundamental theory 
of the birth of the Universe, observational constraints on depar- 
tures from global isotropy are necessary. General Relativity itself 
places no restrictions upon the topology of the Universe, as was 
recognised very early on (e.g., de Sitter 1917); most proposed 
theories of quantum gravity predict topology-change in the early 
Universe which could be visible at large scales today. 

The Einstein field equations relate local properties of the cur- 
vature to the matter content in spacetime. By themselves they do 
not restrict the global properties of the space, allowing a uni- 
verse with a given local geometry to have various global topolo- 
gies. Friedmann-Robertson- Walker (FRW) models of the uni- 
verse observed to have the same average local properties ev- 
erywhere still have freedom to describe quite diff'erent spaces 
at large scales. Perhaps the most remarkable possibility is that a 
vanishing or negative local curvature (Q^k = 1 - ^tot ^ 0) does 
not necessarily mean that our Universe is infinite. Indeed we can 
still be living in a universe of finite volume due to the global 
topological multi-connectivity of space, even if described by the 
flat or hyperbolic FRW solutions. In particular, quantum fluctua- 
tions can produce compact spaces of constant curvature, both flat 
(e.g., Zeldovich & Starobinskii 1984) and curved (e.g., Coule & 
Martin 2000; Linde 2004), within the inflationary scenario. 

The primary CMB anisotropy alone is incapable of con- 
straining curvature due to the well-known geometrical degen- 
eracy which produces identical small-scale fluctuations when 
the recombination sound speed, initial fluctuations, and comov- 
ing distance to the last scattering surface are kept constant 
(e.g.. Bond et al. 1997; Zaldarriaga & Seljak 1997; Stompor 
& Efstathiou 1999). The present results from Planck (Planck 
Collaboration XVI 2013) can therefore place restrictive con- 
straints on the curvature of the Universe only when combined 
with other data: = -K{RqHq)-^ = -0.0010+^;^^!^ at 95 %, 
considering CMB primary anisotropy and lensing from Planck 
(in the natural units with c = 1 we use throughout). This is equiv- 
alent to constraints on the radius of curvature RqHq > 19 for pos- 
itive curvature (K = -hi) and RqHq > 33 for negative curvature 
(K = -1). CMB primary anisotropy alone gives limits on RqHq 
roughly a factor of two less restrictive (and strongly dependent 
on priors). 

Thus, the global nature of the Universe we live in is still an 
open question and studying the observational efl'ects of a possi- 
ble finite universe is one way to address it. With topology not af- 
fecting local mean properties that are found to be well described 
by FRW parameters, its main observational eff'ect is in setting 
boundary conditions on perturbation modes that can be excited 
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and developed into the structure that we observe. Studying struc- 
ture on the last scattering surface is the best-known way to probe 
the global organisation of our Universe and the CMB provides 
the most detailed and best understood dataset for this purpose. 

We can also relax assumptions about the global structure 
of spacetime by allowing anisotropy about each point in the 
Universe. This yields more general solutions to Einstein's field 
equations, leading to the so-called Bianchi cosmologies. For 
small anisotropy, as demanded by current observations, linear 
perturbation about the standard FRW model may be applied. A 
universal shear and rotation induce a characteristic subdominant, 
deterministic signature in the CMB, which is embedded in the 
usual stochastic anisotropics. The deterministic CMB tempera- 
ture fluctuations that result in the homogenous Bianchi models 
were first examined by Collins & Hawking (1973) and Barrow 
et al. (1985), however no dark energy component was included 
as it was not considered plausible at the time. More recently, 
Jaff'e et al. (2006c), and independently Bridges et al. (2007), ex- 
tended these solutions for the open and flat Bianchi Yllh mod- 
els to include cosmologies with dark energy. Although focus 
is given to temperature signatures here, the induced CMB po- 
larization fluctuations that arise in Bianchi models have also 
been derived recently (Pontzen & Challinor 2007; Pontzen 2009; 
Pontzen & Challinor 2011). 

In this paper, we will explicitly consider models of global 
topology and anisotropy. In a chaotic inflation scenario, how- 
ever, our post-inflationary patch might exhibit large-scale local 
topological features ("handles" and "holes") the can mimic a 
global multiply-connected topology in our observable volume. 
Similarly, it might also have residual shear or rotation which 
could mimic the properties of a global Bianchi spacetime. 

Planck's ability to discriminate and remove large-scale astro- 
physical foregrounds (Planck Collaboration XII 2013) reduces 
the systematic error budget associated with measurements of the 
CMB sky significantly. Planck data therefore allow refined lim- 
its on the scale of the topology and the presence of anisotropy. 
Moreover, previous work in this field has been done by a wide 
variety of authors using a wide variety of data (e.g., COBE, 
WMAP 1-year, 3 -year, 5 -year, etc.) and in this work we perform 
a coherent analysis. 

In Sect. 2, we discuss previous attempts to limit the topology 
and global isotropy of the Universe. In Sect. 3 we discuss the sig- 
nals induced in topologically non-trivial and Bianchi universes. 
In Sect. 4 the Planck data we use in the analysis are presented, 
and in Sect. 5 the methods we have developed to detect those 
signals are discussed. We apply those methods in Sect. 6 and 
discuss the results in Sect. 7. 

2. Previous results 

The first searches for non-trivial topology on cosmic scales 
looked for repeated patterns or individual objects in the distri- 
bution of galaxies (Sokolov & Shvartsman 1974; Fang & Sato 
1983; Fagundes & Wichoski 1987; Weatherley et al. 2003). 
The last scattering surface from which the CMB is released 
represents the most distant source of photons in the Universe, 
and hence the largest scales with which we could probe the 
topology of the Universe. This first became possible with the 
DMR instrument on the COBE satellite (Bennett et al. 1996): 
various searches found no evidence for non-trivial topologies 
(e.g., Starobinskij 1993; Sokolov 1993; Stevens et al. 1993; de 
Oliveira-Costa & Smoot 1995; Levin et al. 1998; Bond et al. 
1998, 2000b; Rocha et al. 2004), but sparked the creation of ro- 
bust statistical tools, along with greater care in the enumeration 
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of the possible topologies for a given geometry (see, for exam- 
ple, Lachieze-Rey & Luminet 1995 and Levin 2002 for reviews). 
With data from the WMAP satellite (Jarosik et al. 2011), these 
theoretical and observational tools were applied to a high-quality 
dataset for the first time. Luminet et al. (2003) and Caillerie et al. 
(2007) claimed the low value of the low multipoles (compared to 
standard ACDM cosmology) as evidence for missing large-scale 
power as predicted in a closed universe with a small fundamental 
domain. However, searches in pixel space (Cornish et al. 2004; 
Niarchou et al. 2004; Bielewicz & Riazuelo 2009; Dineen et al. 
2005) and in harmonic space (Kunz et al. 2006) determined that 
this was an unlikely explanation for the low power. Bond et al. 
(1998, 2000a) and Riazuelo et al. (2004a,b) presented some of 
the mathematical formalism for the computation of the correla- 
tions induced by topology in a form suitable for use in cosmo- 
logical calculations. Phillips & Kogut (2006) presented efficient 
algorithms for the computation of the correlation structure of the 
flat torus and applied it via a Bayesian formalism to the WMAP 
data; similar computations for a wider range of geometries were 
performed by Niarchou & Jaffe (2007). 

These calculations used a variety of different vintages of the 
COBE and WMAP data, as well as a variety of different sky cuts 
(including the unmasked internal linear combination (ILC) map, 
not originally intended for cosmological studies). Nonetheless, 
none of the pixel- space calculations which took advantage of the 
full correlation structure induced by the topology found evidence 
for a multiply-connected topology with a fundamental domain 
within or intersecting the last scattering surface. Hence in this 
paper we will attempt to corroborate this earlier work and put 
the calculations on a consistent footing. 

The open and flat Bianchi type VII/^ models have been com- 
pared previously to both the COBE (Bunn et al. 1996; Kogut 
et al. 1997) and WMAP (Jaffe et al. 2005, 2006b) data, albeit 
ignoring dark energy, in order to place limits on the global rota- 
tion and shear of the Universe. A statistically significant corre- 
lation between one of the Bianchi VII/^ models and the WMAP 
ILC map (Bennett et al. 2003) was first detected by Jaffe et al. 
(2005). However, it was noted that the parameters of this model 
are inconsistent with standard constraints. Nevertheless, when 
the WMAP ILC map was "corrected" for the best-fit Bianchi 
template, some of the so-called "anomalies" reported in WMAP 
data disappear (Jaffe et al. 2005, 2006b; Cayon et al. 2006; 
McEwen et al. 2006). A modified template fitting technique was 
performed by Land & Magueijo (2006) and, although a statis- 
tically significant template fit was not reported, the correspond- 
ing "corrected" WMAP data were again free of many large scale 
"anomalies". Due to the consequent renewed interest in Bianchi 
models, solutions to the CMB temperature fluctuations induced 
in Bianchi VII/^ models when incorporating dark energy were 
since derived by Jaffe et al. (2006c) and Bridges et al. (2007). 
Nevertheless, the cosmological parameters of the Bianchi tem- 
plate embedded in WMAP data in this setting remain inconsis- 
tent with constraints from the CMB alone (Jaffe et al. 2006a,c). 
A Bayesian analysis of Bianchi VII/^ models was performed by 
Bridges et al. (2007) using WMAP ILC data to explore the joint 
cosmological and Bianchi parameter space via Markov chain 
Monte Carlo sampling, where it was suggested that the CMB 
"cold spot" (Vielva et al. 2004; Cruz et al. 2006; Vielva 2010) 
could be driving evidence for a Bianchi component (Bridges 
et al. 2008). Recently, this Bayesian analysis has been revisited 
by McEwen et al. (2013) to handle partial-sky observations and 
to use nested sampling methods (Skilling 2004; Feroz & Hobson 
2008; Feroz et al. 2009). McEwen et al. (2013) conclude that 



WMAP 9-year data do not favour Bianchi VII/^ cosmologies over 
ACDM. 



3. CMB correlations in anisotropic and 
multiply-connected universes 

3.1. Topology 

All FRW models can describe multi-connected universes. In the 
case of flat space, there are a finite number of compactifications, 
the simplest of which are those of the torus. All of them have 
continuous parameters that describe the length of periodicity in 
some or all directions (e.g., Riazuelo et al. 2004b). In a space of 
constant non-zero curvature the situation is notably diff'erent — 
the presence of a length scale (the curvature radius 7?o) precludes 
topological compactification at an arbitrary scale. The size of the 
space must now reflect its curvature, linking topological proper- 
ties to Iltot = 1 -^K- In the case of hyperbolic spacetimes, the list 
of possible compact spaces of constant negative curvature is still 
infinite, but discrete (Thurston 1982), while in the positive curva- 
ture spherical space there is only a finite set of well-proportioned 
possibilities (i.e., those with roughly comparable sizes in all di- 
rections; there are also the denumerably infinite lens and prism 
topologies) for a multi-connected space (e.g., Gausmann et al. 
2001; Riazuelo et al. 2004a). 

The eff'ect of topology is equivalent to considering the full 
simply-connected three-dimensional spatial slice of the space- 
time (known as the covering space) as being filled with repeti- 
tions of a shape which is finite in some or all directions {\ht fun- 
damental domain) — by analogy with the two-dimensional case, 
we say that the fundamental domain tiles the covering space. 
For the flat and hyperbolic geometries, there are infinite copies 
of the fundamental domain; for the spherical geometry, with a 
finite volume, there is a finite number of tiles. Physical fields re- 
peat their configuration in every tile, and thus can be viewed as 
defined on the covering space but subject to periodic boundary 
conditions. Topological compactification always break isotropy, 
and for some topologies also the global homogeneity of physical 
fields. Positively curved and flat spaces studied in this paper are 
homogeneous, however hyperbolic multi-connected spaces are 
never homogeneous. 

The primary observable efl'ect of a multi-connected universe 
is the existence of directions in which light could circumnavi- 
gate the space in cosmological time more than once, i.e., the ra- 
dial distance ;^rec to the surface of last scattering exceeds the size 
of the universe. In these cases, the surface of last scattering can 
intersect the (notional) edge of a fundamental domain. At this in- 
tersection, we can view the same spacetime event from multiple 
directions — conversely, it appears in difl'erent directions when 
observed from a single point. 

Thus, temperature perturbations in one direction, T{n), be- 
come correlated with those in another direction T{m) by an 
amount that difl'ers from the usual isotropic correlation function 
C{6), where 6 denotes the angle between h and m. Considering 
a pixelized map, this induces a correlation matrix Cpp> which 
depends on quantities other than the angular distance between 
pixels p and p\ This break from statistical isotropy can there- 
fore be used to constrain topological models. Hence, we need to 
calculate the pixel-space correlation matrix or its equivalent in 
harmonic space. 

In this paper we consider the following topologies: a) 
toroidal flat models with equal-length compactification size L in 
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Table 1: Parameters of analysed curved spaces. 



Spherical Hyperbolic 





Dodecahedral 


Truncated Cube 


Octahedral 


m004(-5,l) 


v3543(2,3) 




0.16 


0.41 


0.82 


0.98 


6.45 




0.31 (;r/10) 


0.39 (;r/8) 


0.45 


0.54 


0.89 




0.37 


0.56 


0.56 


0.64 


1.22 




0.40 


0.58 


0.79 (7T/4) 


0.75 


1.33 



three directions, denoted T[L, L, L];^ b) toroidal flat models with 
different compactification lengths, parametrized by L^^Ly^L^, 
denoted T[Lx,Ly,Lz]; c) three major types of single-action 
positively curved spherical manifolds with dodecahedral, trun- 
cated cubical and octahedral fundamental domains (/*, O*, T* 
compactification groups correspondingly, see Gausmann et al. 
2001); and d) two sample negative curvature hyperbolic spaces, 
m004(-5,l) being one of the smallest known compact hyper- 
bolic spaces as well as the relatively large v3543(2,3).^ Scales of 
fundamental domains of compactified curved spaces are fixed in 
the units of curvature and are summarised in Table 1 , where we 
quote the volume "V, radius of the largest sphere that can be in- 
scribed in the domain "Ri (equal to the distance to the nearest face 
from the origin of the domain), the smallest sphere in which the 
domain can be inscribed "R^ (equal to the distance to the farthest 
vertex), and the intermediate scale that is taken to be the dis- 
tance to the edges for spherical spaces and the "spine" distance 
for hyperbolic topologies. For the cubic torus with edge length L, 
these lengths are = L/2, = V2L/2 and K = V3L/2. The 
ratio "Ru/^ is a good indicator of the shape of the fundamental 
domain. Note that when;^rec is less than Ri, multiple images on 
large scales are not present, although the Cpp' correlation matrix 
is still modified versus the singly-connected limit. The eff'ects of 
topology usually become strong when Xyqc exceeds the interme- 
diate Rjn- 

3.1 .1 . Computing correlation matrices 

The CMB temperature pixel-pixel correlation matrix is a dou- 
ble radial integral of the ensemble average of the product of the 
source functions that describe the transport of photons through 
the universe from the last scattering surface to the observer: 

Cpp'= dx \ dx'{S{xqp)S{x'qp')), (1) 
Jo Jo 

where qp and qp' are unit vectors that point at pixels p and p' 
on the sky, and x and are proper distances along radial rays 
pointing towards the last scattering surface. 

Two techniques have been developed to compute the CMB 
correlation function for multiply-connected universes. In one ap- 
proach, one constructs the orthonormal set of basis functions 
that satisfy the boundary conditions imposed by compactifi- 
cation (eigenfunctions of the Laplacian operator furnish such 
a basis), and assembles the spatial correlation function of the 
source {S(xqp)S(x' ^p')) from such a basis (Cornish & Spergel 

^ In a slight abuse of notation, the lengths Li will be given in units of 
in r[Li , L2, L3], but in physical units elsewhere. 

^ The nomenclature for hyperbolic spaces follows J. Weeks' cen- 
sus, as incorporated in the freely available SnapPea software, http : 
//www. geometrygames.org/SnapPea; see also Thurston & Levy 
(1997). 



1999; Lehoucq et al. 2002). In the other approach, one ap- 
plies the method of images to create the compactified version 

of {S(xqp)^(x' ^p')) froi^ the one computed on the universal 
covering space by resumming the latter over the images of the 
3D spatial positions (Bond et al. 1998, 2000a,b): 

{Sixqp)S{x'qp^)f = Yj<'^(xqp)7[S(7\x' qp'W^ (2) 

where the superscripts c and u refer to the quantity in the 
multiply-connected space and its universal cover, respectively. 
Tilde refers to the need for sum regularization in the models with 
an infinite set of images, e.g., hyperbolic and flat toroidal ones. F 
is the discrete subgroup of motions which defines the multiply- 
connected space and y[x] is the spatial point on the universal 
cover obtained by the action of the motion y g F on the point x. 
The important issue to note is that one needs to implement the 
action of the motion y on the source function unless all the terms 
in the source function are scalar quantities (which is the case if 
one limits consideration to Sachs- Wolfe terms) when the action 
is trivial. 

Both methods are general, but have practical considerations 
to take into account when one increases the pixel resolution. For 
computing Cpp^ up to the resolution corresponding to harmonic 
mode £ ^ 40 both methods have been tested and were found to 
work equally well. In this paper we employ both approaches. 

The main effect of the compactification is that Cpp' is no 
longer a function of the angular separation between the pixels 
p and p' only, due to the lack of global isotropy. In harmonic 
space the two-point correlation function of the CMB is given by 

= {Cl£m<^''^'ni') ^ C^da'^mm' , (3) 

where 6^' is the Kronecker delta symbol and a^m are the spher- 
ical harmonic coefficients of the temperature on the sky when 
decomposed into the spherical harmonics Y^m^q) by 

r(^) = ^a,^y,^($). (4) 

Note that the two-point correlation function C™' is no longer 
diagonal, nor is it m-independent, as in an isotropic universe. 

A flat universe provides an example when the eigenfunctions 
of the Laplacian are readily available in a set of plane waves. 
The topological compactification in the flat space discretizes the 
spectrum of the wavevector magnitudes and selects the subset 
of allowed directions. For example, for a toroidal universe the 
length of the fundamental cell needs to be an integer multiple 
of the wavelength of the modes. We therefore recover a discrete 
sum over modes A:n = (27r/L)n for n = (nx,ny,nz) a triplet of 
integers, instead of an integral over k, 

C^/ oc r d^kAi(k,AT])Ai^(k,AT])P(k) 
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Fig. 1: The top row shows the correlation structure (i.e., a sin- 
gle row of the correlation matrix) of a simply-connected uni- 
verse with isotropic correlations. For subsequent rows, the left 
and middle column show positively curved multiply-connected 
spaces (left: dedocahedral, middle: octahedral) and the right col- 
umn shows equal sided tori. The upper row of three maps cor- 
responds to the case when the size of the fundamental domain 
is of the size of the diameter to the last scattering surface and 
hence the first evidence for large angle excess correlation ap- 
pears. Subsequent rows correspond to decreasing fundamental 
domain size with respect to the last scattering diameter, with pa- 
rameters roughly chosen to maintain the same ratio between the 
models. 



2^ ^,{kn, Ari)AAkn, A7])P(kn)Y,^(n)Y;^Xn) , 



(5) 



where Af(k, At/) is the radiation transfer function (e.g.. Bond & 
Efstathiou 1987; Seljak 1996). We refer to the cubic torus with 
three equal sides as the T3 topology; it is also possible for the 
fundamental domain to be compact in only two spatial dimen- 
sions (e.g., the so-called T2 "chimney" space) or one (the Tl 
"slab", similar to the "lens" spaces available in manifolds with 
constant positive curvature) in which case the sum is replaced by 
an integral in those directions. These models serve as approxi- 
mations to modifications to the local topology of the global man- 
ifold (albeit on cosmological scales): for example, the chimney 
space can mimic a "handle" connecting diff'erent regions of an 
approximately flat manifold. 

In Fig. 1 we show rows of the pixel-space correlation matrix 
for a number of multiply-connected topologies as a map, show- 
ing the magnitude of the correlation within a particular pixel. 
For the simply-connected case, the map simply shows the same 
information as the correlation function C(0); for the topologi- 
cally non-trivial cases, we see the correlations depend on dis- 
tance and direction and difl'er from pixel to pixel (i.e., from row 
to row of the matrix). In Fig. 2 we show example maps of CMB 
anisotropics in universes with these topologies, created by direct 
realisations of Gaussian fields with the correlation matrices of 
Fig. 1. 

3.2. Bianchi 

Bianchi cosmologies include the class of homogeneous but 
anisotropic cosmologies, where the assumption of isotropy about 
each point in the Universe is relaxed. For small anisotropy, as 



Fig. 2: Random realisations of temperature maps for the models 
in Fig. 1. The maps are smoothed with a Gaussian filter with 
fufl- width-half-maximum FWHM = 640'. 



demanded by current observations, linear perturbation about the 
standard FRW model may be applied, leading to a subdominant, 
deterministic contribution to the CMB fluctuations. In this set- 
ting CMB fluctuations may be viewed as the sum of a determin- 
istic Bianchi contribution and the usual stochastic contribution 
that arises in the ACDM model. The deterministic CMB temper- 
ature fluctuations that result in the Bianchi models were derived 
by Barrow et al. (1985), although no dark energy component was 
included. More recently, Jafl'e et al. (2006c), and independently 
Bridges et al. (2007), extended these solutions for the open and 
flat Bianchi VII/^ models to include cosmologies with dark en- 
ergy. We defer the details of the CMB temperature fluctuations 
induced in Bianchi models to these works and give only a brief 
description here. 

Bianchi Yllh models describe a universe with overall ro- 
tation, parameterized by an angular velocity, oj, and a three- 
dimensional rate of shear, parameterized by the antisymmetric 
tensor ct/^; we take these to be relative to the z axis. The model 
has a free parameter, first identified by Collins & Hawking 
(1973), describing the comoving length-scale over which the 
principal axes of shear and rotation change orientation. The ra- 
tio of this length scale to the present Hubble radius is typically 
denoted x, which defines the h parameter of type Yllh models 
through (Barrow et al. 1985) 



X = 



h 



(6) 



where the total energy density Iltot = + ^a- The parameter 
X acts to change the "tightness" of the spiral-type CMB tem- 
perature contributions that arise due to the geodesic focusing of 
Bianchi VII/^ cosmologies. The shear modes CTij of combinations 
of orthogonal coordinate axes are also required to describe a 
Bianchi cosmology. The present dimensionless vorticity (oj/H)o 
may be related to the dimensionless shear modes (crij/H)o by 
(Barrow et al. 1985) 



a- 



(1 -h /z)i/2(i + 9/i) 

6x^Qfa 



1/2 



(7) 



where H is the Hubble parameter. Throughout we assume equal- 
ity of shear modes cr = an = (Xn (cf. Jaff'e et al. 2005). The 
amplitude of the deterministic CMB temperature fluctuations in- 
duced in Bianchi Yllh cosmologies may be characterised by ei- 



5 



Planck Collaboration: Planck 2013 results. XXVL Background geometry and topology of the Universe 






Fig. 3: Simulated deterministic CMB temperature contribu- 
tions in Bianchi Nllh cosmologies for varying x and Qtot 
(left-to-right Iltot e {0.10,0.30,0.95}; top-to-bottom x e 
{0.1,0.3,0.7, 1.5,6.0}). In these maps the swirl pattern typical 
of Bianchi-induced temperature fluctuations is rotated from the 
South pole to the Galactic centre for illustrational purposes. 



ther (cr/H)o or (oj/H)o since these parameters influence the am- 
plitude of the induced temperature contribution only and not its 
morphology. The handedness of the coordinate system is also 
free in Bianchi VII/^ models, hence both left- and right-handed 
models arise. Since the Bianchi-induced temperature fluctua- 
tions are anisotropic on the sky the orientation of the result- 
ing map may vary also, introducing three additional degrees-of- 
freedom. The orientation of the map is described by the Euler 
angles^ (a,p,y), where for (a,/3,y) = (0°,0°,0°) the swirl pat- 
tern typical of Bianchi templates is centred on the South pole. 

Examples of simulated Bianchi VII/j CMB temperature maps 
are illustrated in Fig. 3 for a range of parameters. In the anal- 
ysis performed herein the BIANCHI 2^ (McEwen et al. 2013) 
code is used to simulate the temperature fluctuations induced 
in Bianchi VII/^ models. Bianchi VII/^ models induce only large 
scale temperature fluctuations in the CMB and consequently 
Bianchi maps have a particularly low band-limit, both globally 
and azimuthally (i.e., in both £ and m in spherical harmonic 
space; indeed, only those harmonic coefficients with m = ±1 
are non-zero). 



4. Data description 

We use Planck maps that have been processed by the 
various component-separation pipelines described in Planck 
Collaboration XII (2013). The methods produce largely consis- 
tent maps of the sky, with detailed diff'erences in pixel intensity, 
noise properties, and masks. Here, we consider maps produced 
by the Commander -Ruler, NILC, SMICA and SEVEM methods. 
Each provides its own mask and we also consider the conserva- 
tive common mask. 

We note that because our methods rely on rather intensive 
pixel- or harmonic- space calculations, in particular considering 




^ The active zyz Euler convention is adopted, corresponding to the 
rotation of a physical body in a fixed coordinate system about the z, y 
and z axes by y, /3 and a respectively. 

^ http://www.jasonmcewen.org/ 



Fig. 4: The mask (/sky = 0.76) used in the matched circles anal- 
ysis. 



a full set of three-dimensional orientations and, for the likeli- 
hood methods, manipulation of an anisotropic correlation ma- 
trix, computational efficiency requires the use of data degraded 
from the native HEALPix (Gorski et al. 2005) A^side = 2048 
resolution of the Planck maps. Because the signatures of ei- 
ther a multiply-connected topology or a Bianchi model are most 
prominent on large angular scales, this does not result in a sig- 
nificant loss of ability to detect and discriminate amongst the 
models (see Sect. 5.3). 

The topology analyses both rely on degraded maps and 
masks. The matched-circles method smooths with a 30' 
Gaussian filter and degrades the maps to A^side = 512, and uses 
a mask derived from the SEVEM component separation method 
(Fig. 4). Because the performance of the matched-circles statistic 
can be significantly degraded by the point source cut, we mask 
only those point sources from the full-resolution /sky = 0.73 
SEVEM mask with amplitude, after smoothing and extrapolation 
to the 143 or 217 GHz channels, greater than the faintest source 
originally detected at those frequencies. The mask derived in this 
way retains /ky = 0.76 of the sky. 

The likelihood method smooths the maps and masks with an 
11° Gaussian filter and then degrades them to A^side = 16 and 
conservatively masks out any pixel with more than 10 % of its 
original subpixels masked. At full resolution, the common mask 
retains a fraction /ky = 0.73 of the sky, and /ky = 0.78 when 
degraded to A^side = 16 (the high-resolution point-source masks 
are largely filled in the degraded masks). The Bianchi analysis 
is performed in harmonic space, and so does not require explicit 
degradation in pixel space. Rather, the data are transformed at 
full resolution into harmonic space and considered only up to a 
specified maximum harmonic ^, where correlations due to the 
mask are taken into account. 

Diff'erent combinations of these maps and masks are used 
to discriminate between the topological and anisotropic models 
described in Sect. 3. 



5. Methods 

5.1. Topology: circles in the sky 

The first set of methods, exemplified by the circles-in-the-sky 
of Cornish et al. (1998), involves a frequentist analysis using a 
statistic which is expected to diff'er between the models exam- 
ined. For the circles, this uses the fact that the intersection of the 
topological fundamental domain with the surface of last scat- 
tering is a circle, which one potentially views from two diff'er- 
ent directions in a multiply-connected universe. Of course, the 
matches are not exact due to noise, foregrounds, the integrated 
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Sachs- Wolfe (ISW) and Doppler effects along the different lines 
of sight. 

By creating a statistic based on the matching of differ- 
ent such circles, we can compare Monte Carlo simulations of 
both a simply-connected, isotropic null model with specific 
anisotropic or topological models. We may then calibrate de- 
tections and non-detections using Monte Carlo simulations. In 
principle, these simulations should take into account the com- 
plications of noise, foreground contributions, systematics, the 
ISW and Doppler effects. However, they do not include gravi- 
tational lensing of the CMB as the lensing deflection angle is 
small compared to the minimal angular scale taken into account 
in our analysis. Note that the null test is generic (i.e., not tied 
to a specific topology) but any detection must be calibrated with 
specific simulations for a chosen topology or anisotropic model. 
A very similar technique can be used for polarisation by taking 
into account the fact that the polarisation pattern itself is now 
not directly repeated, but rather that the underlying quadrupole 
radiation field around each point on the sky is now seen from 
different directions (Bielewicz et al. 2012). These methods have 
been applied successfully to COBE DMR and WMAP data, and 
have recently been shown to be feasible for application to Planck 
data (Bielewicz et al. 2012). 

If light had sufficient time to cross the fundamental domain, 
an observer would see multiple copies of a single astronomical 
object. To have the best chance of seeing "around the Universe" 
we should look for multiple images of the furthest reaches of the 
Universe. Searching for multiple images of the last scattering 
surface — the edge of the visible Universe — is then a powerful 
way to constrain topology. Because the surface of last scattering 
is a sphere centred on the observer, one can imagine each copy 
of the observer will come with a copy of the last scattering sur- 
face, and if the copies are separated by a distance less than the 
diameter of the last scattering surface, then they will intersect 
along circles. These are visible to both copies of the observer, but 
from opposite sides. The two copies are really one observer, so if 
space is sufficiently small, the CMB radiation from the last scat- 
tering surface will demonstrate a pattern of hot and cold spots 
that matches around the circles. 

The idea of using these circles to study topology is due to 
Cornish et al. (1998). In that work, a statistical tool was devel- 
oped to detect correlated circles in all sky maps of the CMB 
anisotropy — the circle comparison statistic. In our studies we 
will use version of this statistic optimised for the small-scale 
anisotropics as defined by Cornish et al. (2004): 

Stjia, cf>.) = '^^^^-^ , (8) 

i:.N(|Ar,„p + |Ar,;„p) 

where ATi^m and ATj^^ denote the Fourier coefficients of the tem- 
perature fluctuations around two circles of angular radius a cen- 
tered at different points on the sky, / and j, respectively, with rel- 
ative phase 0*. The mth harmonic of the temperature anisotropics 
around the circle is weighted by the factor |m|, taking into ac- 
count the number of degrees of freedom per mode. Such weight- 
ing enhances the contribution of small-scale structure relative 
to large-scale fluctuations and is especially important since the 
large-scale fluctuations are dominated by the ISW effect. This 
can obscure the image of the last scattering surface and reduce 
the ability to recognise possible matched patterns on it. 

The above S ^ statistic corresponds to pair of circles with the 
points ordered in a clockwise direction (phased). For alternative 
ordering, when along one of the circles the points are ordered 



in an anti-clockwise direction (anti-phased), the Fourier coeffi- 
cients ATi^tn are complex conjugated, defining the S~ statistic. 
This allows the detection of both orientable and non-orientable 
topologies. For orientable topologies the matched circles have 
anti-phased correlations while for non-orientable topologies they 
have a mixture of anti-phased and phased correlations. 

The statistic has a range over the interval [-1,1]. Circles that 
are perfectly matched have S = 1, while uncorrelated circles will 
have a mean value ofS =0. Although the statistic can also take 
negative values for the temperature anisotropy generated by the 
Doppler term (Bielewicz et al. 2012), anticorrelated circles are 
not expected for the total temperature anisotropy considered in 
this work. To find matched circles for each radius a, the maxi- 
mum value S^^^(a) = max/j;^^ ^fj(^^ 0*) determined. 

Because general searches for matched circles are computa- 
tionally very intensive, we restrict our analysis to a search for 
pairs of circles centered around antipodal points, so called back- 
to-back circles. As described above, the maps were also down- 
graded to /"/side = 512, which greatly speeds up the computations 
required, but with no significant loss of discriminatory power, as 
seen in Sect. 5.3.1. More details on the numerical implementa- 
tion of the algorithm can be found in the paper by Bielewicz & 
Banday (2011). 

The constraints we will derive concern only topologies that 
predict matching pairs of back- to-back circles. Thus, we cannot 
rule out inhomogeneous spaces for which the relative position of 
the circles depends on the position of the observer in the funda- 
mental polyhedron (Bond et al. 2000b; Riazuelo et al. 2004b). 
The strongest constraints are imposed on topologies predicting 
back- to-back circles in all directions i.e., all the single action 
manifolds, among them tori of any shape and the three spherical 
cases considered in this paper. Weaker constraints are imposed 
on topologies with all back-to-back circles centred on a great cir- 
cle of the celestial sphere such as half-turn, quarter-turn, third- 
turn and sixth-turn spaces, as well as Klein and chimney spaces. 
The statistic can also constrain the multi-connected spaces pre- 
dicting one pair of antipodal matching circles such as Klein or 
chimney spaces with horizontal flip, vertical flip or half-turn and 
slab space translated without screw motion. Other topologies 
catalogued in Riazuelo et al. (2004b) are not constrained by this 
analysis: the Hantzsche-Wendt space; the chimney space with 
half-turn and flip; the generic slab space; the slab space with flip; 
spherical manifolds with double and linked action; and all the 
hyperbolic topologies. Note that we explicitly search for some 
of these cases with the likelihood method discussed in Sect. 5.2 
below. 

Furthermore, even for the topologies that predict matching 
pairs of back-to-back circles, the constraints do not apply to 
those universes for which the orientation of the matched circles 
is impossible to detect due to partial masking on the sky. Because 
of the larger sky fraction removed by the Planck common mask 
than for WMAP this probability is larger for the analysis of the 
Planck maps. Moreover, the smaller fraction of the sky used in 
the search of matched circles results in a false detection level 
larger with our /sky = 0.76 mask than for the /sky = 0.78 7-year 
KQ85 WMAP mask. As a result we obtain weaker — but more 
conservative — constraints on topology than for similar analyses 
of WMAP data (Bielewicz & Banday 201 1). 

To draw any conclusions from an analysis based on the statis- 
tic S^^^(a), it is very important to correctly estimate the thresh- 
old for a statistically significant match of circle pairs. We used 
300 Monte Carlo simulations of CMB maps, described in detail 
in Section 5.3.1, to establish the threshold such that fewer than 
1 % of simulations would yield a false event. 
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5.2. Bayesian analyses 

The second set of methods take advantage of the fact that the 
underlying small-scale physics is unchanged in both anisotropic 
and topological models compared to the standard cosmology, 
and thus a Gaussian likelihood function will still describe the 
statistics of temperature and polarization on the sky, albeit no 
longer with isotropic correlations. When considering specific 
topologies, these likelihood methods instead calculate the pixel- 
pixel correlation matrix. This has been done for various torus 
topologies (which are a continuous family of possibilities) in 
the flat Universe as well as for locally hyperbolic and spher- 
ical geometries (which have a discrete set of possibilities for 
a given value of the curvature). More general likelihood-based 
techniques have been developed for generic mild anisotropics 
in the initial power spectrum (Hanson & Lewis 2009), which 
may have extension to other models. For the Bianchi setting, 
an isotropic zero-mean Gaussian likelihood is recovered by sub- 
tracting a deterministic Bianchi component from the data, where 
the cosmological covariance matrix remains diagonal in har- 
monic space but masking introduces non-diagonal structure that 
must be taken into account. 

Because these methods use the likelihood function directly, 
they can take advantage of any detailed noise correlation infor- 
mation that is available, including any correlations induced by 
the foreground-removal process. We denote the data by the vec- 
tor d, which may be in the form of harmonic coefficients d^m or 
pixel temperatures dp or, in general, coefficients of the tempera- 
ture expansion in any set of basis functions. We denote the model 
under examination by the discrete parameter M, which can take 
on the appropriate value denoting the usual isotropic case, or 
the Bianchi case, or one of the possible multiply-connected uni- 
verses. The continuous parameters of model M are given by 
the vector 0, which for this case we can partition into ©c for 
the cosmological parameters shared with the usual isotropic 
and simply-connected case, and ©a which denotes the param- 
eters for the appropriate anisotropic case, be it a topologically 
non-trivial universe or a Bianchi model. Note that all of the 
anisotropic cases contain "nuisance parameters" which give the 
orientation of either the fundamental domain or the Bianchi tem- 
plate which we can marginalize over as appropriate. 

Given this notation, the posterior distribution for the param- 
eters of a particular model, M, is given by B ayes' theorem: 



Pi©\d,M) 



p(e\M)P(d\e, M) 

P{d\M) 



(9) 



Here, P(0|M) = P(0c,0aI^) is the joint prior probability of 
the standard cosmological parameters ©c and those describing 
the anisotropic universe ©a, P{d\&, M) = X is the likelihood, 
and the normalizing constant P{d\M) is the Bayesian evidence, 
which can be used to compare the models to one another. 

We will usually take the priors to be simple "non- 
informative" distributions (e.g., uniform over the sphere for ori- 
entations, uniform in length for topology scales, etc.) as appro- 
priate. The form of the likelihood function will depend on the 
anisotropic model: for multiply-connected models, the topol- 
ogy induces anisotropic correlations, whereas for the Bianchi 
model, there is a deterministic template, which depends on the 
Bianchi parameters, in addition to the standard isotropic cos- 
mological perturbations. We will assume that any other non- 
Gaussian signal (either from noise or cosmology) is negligible 
(Planck Collaboration XXIII 2013; Planck Collaboration XXIV 
2013) and use an appropriate multivariate Gaussian likelihood. 



Given the signal and noise correlations, and a possible 
Bianchi template, the procedure is similar to that used in stan- 
dard cosmological-parameter estimation, with a few complica- 
tions. Firstly, the evaluation of the likelihood function is compu- 
tationally expensive and usually limited to large angular scales. 
This means that in practice the eff'ect of the topology on the like- 
lihood is usually only calculated on those large scales. Secondly, 
the orientation of the fundamental domain or Bianchi template 
requires searching (or marginalizing) over three additional pa- 
rameters, the Euler angles. 

5.2.1. Topology 

In topological studies, the parameters of the model consist of 
©c, the set of cosmological parameters for the fiducial best- 
fit flat cosmological model, and ©t, the topological parameters 
which include the set of compactification lengths Lx,Ly, for 
flat toroidal model or the curvature parameter ^Ik for curved 
spaces, and a choice of compactification T. In our studies we 
keep ©c fixed, and vary ©t for a select choice of compactifica- 
tions listed in Sect. 3.1. These parameters define the predicted 
two-point signal correlation matrix Cpp^ for each model, which 
are precomputed. Additional internal parameters, including the 
amplitude of the signal A and the angles of orientation of the 
fundamental domain of the compact space relative to the sky (f 
(e.g., parameterized by a vector of the three Euler angles), are 
maximized and/or marginalized over during likelihood evalua- 
tion. 

The likelihood, i.e., the probability to find a temperature data 
map d with associated noise matrix N given a certain topological 
model is then given by 

P(^|C[©c,©T,r],A,^) 

oc ^ exp\--d'(AC-^Mr^d\ . (10) 
VIAC + N| "^1 2 ' 'J 

Working with a cut-sky, it is often easier to start the anal- 
ysis with data and a correlation matrix given in pixel space. 
However, especially in the realistic case of negligible noise on 
large scales, the matrix C -h N is poorly conditioned in pixel 
space, and pixel space evaluation of the likelihood is, as a rule, 
not robust. Indeed, there are typically more pixels than indepen- 
dent modes that carry information about the signal (e.g., even 
in the standard isotropic case, sub-arcminute pixels would not 
be useful due to beam- smoothing; with anisotropic correlations 
and masked regions of the sky, more complicated linear combi- 
nations of pixels even on large scales may have very little signal 
content). Therefore in general we expand the temperature map 
dp, the theoretical correlation matrix Cpp^ and the noise covari- 
ance matrix Npp' in a discrete set of mode functions i/^nip), or- 
thonormal over the pixelized sphere, possibly with weights w(p), 
Z/? ^(p)^n(p)^l'(p) = ^nn', obtaining the coefficients of expan- 
sion 

dn = Y^dp\l/l{p)w{p)\ 

p 

Cnn' = ^Y^Cpp'l//n(p)l//l^(p)w(p)w(p); 

P P' 

Nnn' = ^^C^p'<A«(/^¥:'(/^0w(/7)w(/70. (11) 



Next we select Nj^ such modes for comparison and consider the 
likelihood marginalized over the remainder of the modes 

p(d\C[ec,&T,Tl(f,A)oc 
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VIAC + N|m ( 2^ J 

where C and N are restricted to the A^^ x A^m block of chosen 
modes. Flexibility in choosing mode functions and their number 
A^m is used to achieve the compromise between robust invert- 
ibility of the C + N matrix on the one hand, and the amount of 
discriminating information retained in the data on the other. The 
weights w{p) can be used to improve the accuracy of transforms 
on a pixelized sky. 

For full-sky analysis the natural choice of the mode functions 
is the set of ordinary spherical harmonics Yim{p) which leads to 
standard harmonic analysis 

P(J,^|C[0c,0T,7^],^,A) 

cx— L=exp|-i ^ 4yC + N)-J^,,^,J,,^A. (13) 
VIC + N| [ 2^^^, j 

Here, where we focus on masked data, we have made a some- 
what different choice. As a mode set for comparison we use the 
A^m = 837 largest eigenvectors of the Cpp' matrix, restricted to 
the masked sky, for the fiducial flat isotropic model with best-fit 
parameters ©c- We emphasize that the correlation matrix com- 
puted for this reduced dataset has fewer modes, but contains no 
additional assumptions beyond those of the original Cpp' . 

Since computation of Cpp' matrices for a range of topolog- 
ical models is expensive, we do not aim to determine the full 
Bayesian evidence P(d\T) which would require marginalization 
over all parameters ©c, ©t, amplitude A, and orientation (Euler 
angles) (f, and would in addition be sensitive to the prior proba- 
bilities assumed for the size of the fundamental domain. Instead 
we directly compare the likelihood along the changing set of ©t 
that has as its limit the flat fiducial model defined by ©c- In case 
of toroidal topology such a limit is achieved by taking compact- 
ification lengths to infinity, while for curved models we vary 
in comparison to the flat limit = 0. In the latter case, for 
the spherical spaces we change Qa and Hq together with to 
track the CMB geometrical degeneracy line in which the recom- 
bination sound speed, initial fluctuations, and comoving distance 
to the last scattering surface are kept constant (e.g.. Bond et al. 
1997; Zaldarriaga & Seljak 1997; Stompor & Efstathiou 1999), 
and for hyperbolic spaces we vary while keeping Hq and 
fix - fixed to fiducial values. Note that hyperbolic multi- 
connected spaces, in contrast to tori and the single-action pos- 
itive curvature manifolds considered in this paper, are not only 
anisotropic but also inhomogeneous. Therefore, the likelihood is 
expected to be dependent on the position of the observer. We do 
not study this dependence here. 

For each parameter choice, we find the likelihood at the 
best orientation (f of the topology with respect to the sky after 
marginalizing over the amplitude A of the signal (hence, this can 
be considered a profile likelihood with respect to the orientation 
parameters). This likelihood is compared both with the fiducial 
model applied to the observed temperature map and with the 
likelihood of the topological model applied to the simulated re- 
alization of the isotropic map drawn from the fiducial model. 
Such a strategy is optimized for the detection of topological sig- 
natures. For non-detections, the marginalized likelihood can be a 
better probe of the overall power of the data to reject a non-trivial 
topology, and so for real data below, we also show the likelihood 
marginalized over the orientations (f. We estimate the marginal- 
ized likelihood from the random sample of 10,000 orientations, 
drawn statistically uniformly on the sphere of unit quater- 



nions representing rotations of the fundamental domain relative 
to the observed sky. 

5.2.2. Bianchi 

For the Bianchi analysis the posterior distribution of the parame- 
ters of model M is given by B ayes' Theorem, specified in Eq. 9, 
similar to the topological setting. The approach of McEwen et al. 
(2013) is followed, where the likelihood is made explicit in the 
context of fitting a deterministic Bianchi template embedded in 
a stochastic CMB background, defined by the power spectrum 
Q(©c) for a given cosmological model with parameters ©c- 
The Bianchi VII/^ parameters are denoted ©b . The corresponding 
likelihood is given by 

P(^|©B, ©c) oc I cxp[-x\ec, ©b)/2] , (14) 
V|X(©c)| 

where 

x\ec, ©b) = [d- b(eB)fx-\ec)[d - z^(©b)] (15) 

and d = [dim] and Z>(©b) = {^^m(©B)} are the spherical har- 
monic coefficients of the data and Bianchi template, respectively, 
considered up to the harmonic band-limit ^max- A band-limit of 
^max = 32 is considered in the subsequent analysis since this is 
sufficient to capture the structure of the CMB temperature fluctu- 
ations induced in Bianchi VII/^ models (see, e.g., McEwen et al. 
2006). The likelihood is computed in harmonic space where 
rotations of the Bianchi template can be performed efficiently. 
The covariance matrix X(©c) depends on whether the full- sky 
or partial-sky masked setting is considered. In the full-sky set- 
ting X(©c) = C(©c) as first considered by Bridges et al. (2007), 
where C(©c) is the diagonal CMB covariance matrix with en- 
tries Q(©c) on the diagonal. In the case of a zero Bianchi com- 
ponent, Eq. 14 then reduces to the likelihood function used com- 
monly to compute parameter estimates from the power spectrum 
estimated from CMB data (e.g. Verde et al. 2003). In the masked 
setting considered subsequently, X(©c) = C(©c) + M, where M 
is the non-diagonal mask covariance matrix, as considered by 
McEwen et al. (2013). The;^^ of the likelihood for the Bianchi 
case diff'ers from the topology case by the nonzero Bianchi tem- 
plate b and the use of a correlation matrix M to account for the 
presence of the mask. 

In the most physically motivated scenario, the Bianchi and 
cosmological parameters are coupled (e.g. the total density of 
the Bianchi and standard cosmological model are identical). 
However, it is also interesting to consider Bianchi templates as 
phenomenological models with parameters decoupled from the 
standard cosmological parameters, particularly for comparison 
with previous studies. Both scenarios are considered in the sub- 
sequent analysis. In the decoupled scenario a flat cosmologi- 
cal model is considered, whereas in the decoupled scenario an 
open cosmological model is considered to be consistent with the 
Bianchi VII/^ model; we label these models the flat-decoupled- 
Bianchi model and the the open-coupled-Bianchi model, respec- 
tively. 

To determine whether the inclusion of a Bianchi component 
better describes the data the Bayesian evidence is examined, as 
given by 

E = P(d\M) = J deP(d\e,M)P(e\M) . (16) 

Using the Bayesian evidence to distinguish between models nat- 
urally incorporates Occam's razor, trading ofl' model simplicity 
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and accuracy. In the absence of any prior information on the pre- 
ferred model, the Bayes factor given by the ratio of Bayesian 
evidences (i.e., E1/E2) is identical to the ratio of the model prob- 
abilities given the data. The Bayes factor is thus used to distin- 
guish models. The Jeffreys scale (Jeffreys 1961) is often used as 
a rule-of-thumb when comparing models via their Bayes factor. 
The log-Bayes factor Aln£^ = ln(Ei/E2) (also called the log- 
evidence difference) represents the degree by which the model 
corresponding to Ei is favoured over the model correspond- 
ing to £"2, where: < Aln£^ < 1 is regarded as inconclusive; 
1 < AlnE < 2.5 as significant; 2.5 < AlnE < 5 as strong; 
and Aln£^ > 5 as conclusive (without loss of generality we have 
assumed Ei > E2). For reference, a log-Bayes factor of 2.5 cor- 
responds to odds of 1 in 12, approximately, while a factor of 5 
corresponds to odds of 1 in 150, approximately. 

The ANICOSMO^ code (McEwen et al. 2013) is used to per- 
form a Bayesian analysis of Bianchi VII/^ models, which in turn 
uses the public MultiNest^ code (Feroz & Hobson 2008; Feroz 
et al. 2009) to sample the posterior distribution and compute evi- 
dence values by nested sampling (Skilling 2004). We sample the 
parameters describing the Bianchi VII/^ model and those describ- 
ing the standard cosmology simultaneously. 

5.3. Simulations and Validation 
5.3.1. Topology 

Circles-in-the-Sky Before beginning the search for pairs of 
matched circles in the Planck data, we validate our algorithm 
using simulations of the CMB sky for a universe with 3 -torus 
topology for which the dimension of the cubic fundamental do- 
main is L = 2//q\ and with cosmological parameters corre- 
sponding to the AC DM model (see Komatsu et al. 2011, Table 
1) determined from the 7-year WMAP results combined with the 
measurements of the distance from the baryon acoustic oscilla- 
tions and the Hubble constant. We performed simulations com- 
puting directly the a^m coefficients up to the multipole of order 
^max = 500 as described in Bielewicz & Banday (201 1) and con- 
volving them with the same smoothing beam profile as used for 
the data, i.e., a Gaussian beam with 30' FWHM. In particular, 
we verified that our code is able to find all pairs of matched 
circles in such a map. The map with marked pairs of matched 
circles with radius a ^ 24° and the statistic S:^^^(a) for the map 
are shown in Fig. 5 and Fig. 6, respectively. Note that the peak 
amplitudes in the statistic, corresponding to the temperature cor- 
relation for matched circles, decrease with radius of the circles. 
Cornish et al. (2004) noted that this is primarily caused by the 
Doppler term, which becomes increasingly anticorrelated for cir- 
cles with radius smaller than 45°. 

The intersection of the peaks in the matching statistic with 
the false detection level estimated for the CMB map correspond- 
ing to the simply-connected universe defines the minimum ra- 
dius of the correlated circles which can be detected for this map. 
The height of the peak with the smallest radius seen in Fig. 6 
indicates that the minimum radius is about aj^in ~ 20° . 

For the Monte Carlo simulations of the CMB maps for the 
simply-connected universe we used the same cosmological pa- 
rameters as for the multi-connected universe, i.e., corresponding 
to the ACDM model determined from the 7-year WMAP results. 
The maps were also convolved with the same beam profile as for 
the simulated map for the 3 -torus universe and data, as well as 



^ http://www.jasonmcewen.org/ 

^ http : //www . mrao . cam . ac . uk/sof tware/multinest/ 




Fig. 5: A simulated map of the CMB sky in a universe with a 
T[2, 2, 2] toroidal topology. The dark circles show the locations 
of the same slice through the last scattering surface seen on op- 
posite sides of the sky. They correspond to matched circles with 
radius a ^ 24° . 
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Fig. 6: An example of the S^^^ statistic as a function of circle ra- 
dius a for a simulated CMB map (shown in Fig. 5) of a universe 
with the topology of a cubic 3-torus with dimensions L - 2H~^ 
(solid line). The dash-dotted line show the false detection level 
established such that fewer than 1 % out of 300 Monte Carlo 
simulations of the CMB map, smoothed and masked in the same 
way as the data, would yield a false event. 



masked with the same cut (Fig. 4) used for the analysis of data. 
The false detection threshold was established such that fewer 
than 1 % of 300 Monte Carlo simulations would yield a false 
event. 



Bayesian Analysis Because of the expense of the calculation 
of the correlation matrix, we wish to limit the number of three- 
dimensional wavevectors k we consider, as well as the number 
of spherical harmonic modes and finally the number of dif- 
ferent correlation matrices as a whole. We need to ensure that 
the full set of matrices C™' that we calculate contains all of the 
available information on the correlations induced by the topol- 
ogy in a sufficiently fine-grained grid. For this purpose, we con- 
sider the Kullback-Leibler (KL) divergence as a diagnostic (see, 
e.g., Kunz et al. 2006, 2008, for applications of the KL diver- 
gence to topology). The KL divergence between two probability 
distributions pi(x) and P2(x) is given by 



= I PiW^^ — —ax. 



Piix) 



(17) 
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If the two distributions are Gaussian with correlation matrices 
Ci and C2, this expression simplifies to 



(18) 



and is thus a measure of the discrepancy between the correlation 
matrices. The KL divergence can be interpreted as the ensemble 
average of the log-likelihood-ratio AlnX between realizations of 
the two distributions. Hence, they enable us to probe the ability 
to tell if, on average, we can distinguish realizations of pi from 
a fixed p2 without having to perform a brute-force Monte Carlo 
integration. Thus, the KL divergence is related to ensemble aver- 
ages of the likelihood-ratio plots that we present for simulations 
(Fig. 11) and real data (Sect. 6), but does not depend on simu- 
lated or real data. 

We first use the KL divergence to determine the size of the 
fundamental domain which we can consider to be equivalent to 
the simply-connected case (i.e., the limit in which all dimen- 
sions of the fundamental domain go to infinity). We note that 
in our standard ACDM model, the distance to the surface of 
last scattering is Xrtc ~ 3.1416(//o)~^- We would naively ex- 
pect that as long as the sphere enclosing the last scattering sur- 
face can be enclosed by the fundamental domain (L = ^Xvtc), 
we would no longer see the eff'ects of non-trivial topology. 
However, because the correlation matrix includes the full three- 
dimensional correlation information (not merely the purely ge- 
ometrical eff'ects of completely correlated points) we would see 
some long- scale correlation eff'ects even for larger fundamental 
domains. In Fig. 7 we show the KL divergence (as a function 
of (LHo)~^ so that the simply-connected limit L ^ 00 is at a fi- 
nite position) for the T[L, L, L] (cubic), T[L, L, 7] (chimney) and 
r[L, 7,7] (slab) spaces and show that it begins to level off' for 
(LHo)~^ < 1/5, although these topologies are still distinguish- 
able from the r [7, 7, 7] torus which is yet closer to the value for 
a simply-connected universe Jkl[7, 7, 7] ^ 1.1. These figures, as 
well as the likelihoods computed on simulations and data, show 
steps and other structures on a variety of scales corresponding 
to the crossing of the diff'erent length scales of the fundamen- 
tal domain "Rj^, and "Ri crossing the last scattering surface; 
smaller fundamental domains with longer intersections with the 
last scattering surface are easier to detect. 

Computational limitations further prevent us from calculat- 
ing the likelihood at arbitrary values of the fundamental do- 
main size parameters. We must therefore ensure that our coarse- 
grained correlation matrices are sufficient to detect a topology 
even if it lies between our gridpoints. In Fig. 8 we show the 
KL divergence as a function of the size of the fundamental 
domain, relative to various models, both aligned with our grid 
(LHq = 4.5) and in between our grid points (LHq = 5.25). We 
see that the peak is wide enough that we can detect a peak within 
SLHq ~ 0. 1 of the correct value. We also show that we can detect 
anisotropic fundamental domains even when scanning through 
cubic tori: we show a case which approximates a "chimney" uni- 
verse with one direction much larger than the distance to the last 
scattering surface. 

Because our topological analyses do not simultaneously vary 
the background cosmological parameters along with those de- 
scribing the topology, we also probe the sensitivity to the cos- 
mology. In Fig. 9 we show the eff'ect of varying the ffducial cos- 
mology from the Planck Collaboration XVI (2013) best-fft val- 
ues to those reported by WMAP (Komatsu et al. 2011).^ We see 




Fig. 7: The KL divergence computed for torus models as a func- 
tion of the (inverse) length of a side of the cube. r[Li,L2,L3] 
refers to a torus with edge lengths L/. 
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^ We use the wmap7+bao+h0 results from http://lambda.gsfc. 
nasa. gov. 



Fig. 8: The KL divergence between a supposed correct model 
and other models. We show diff'erences of cubic tori with respect 
to models with (LHq)'^ = 1 /4.5 ^ 0.22 (aligned with our grid of 
models), (LHq)-^ = 1/5.25 ^ 0.19 (in between the gridpoints) 
and and a r[5, 5, 7] chimney model with (LHq)'^ = 1/5 in two 
directions and (LHo)'^ = I/7 ^ 0.14 in the third. 



that this induces a small bias of SLHq ^0.2 but does not hinder 
the ability to detect a non-trivial topology. 

We have also directly validated the topological Bayesian 
techniques with simulations. In Fig. 10 we show the log- 
likelihood for the above T[2, 2, 2] simulations as a function of 
two of the Euler angles, maximized over the third. We ffnd a 
strong peak at the correct orientation, with a multiplicity due 
to the degenerate orientations corresponding to the faces of the 
cube (there are peaks at the North and South poles, which are dif- 
ffcult to see in this projection). Note that the peaks correspond to 
ratios of more than exp(700) compared to the relatively smooth 
minima elsewhere. 

In Fig. 11 we also test the ability of the Bayesian likeli- 
hood technique to detect the compactiffcation of the space in the 



11 



Planck Collaboration: Planck 2013 results. XXVL Background geometry and topology of the Universe 




I 0.15 0.20 0.25 0.30 



Fig. 9: The KL divergence between a model generated with the 
WMAP best-fit cosmological parameters as a background cos- 
mology and a r[5,5,5] cubic torus topology with respect to a 
Planck best-fit cosmology and a varying cubic topology. 
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Fig. 10: The log-likelihood with respect to the peak as a function 
of the orientation of the fundamental r[2, 2, 2] torus domain for 
the simulations. The third Euler angle is marginalized over. We 
see peaks at the orientations corresponding to the six faces of 
the cubic fundamental domain (there are peaks at the North and 
South poles, which are difl&cult to see in this projection). 



simulated temperature realizations drawn from the dodecahedral 
closed model. For curved geometries, the size of the fundamen- 
tal domain is fixed with respect to the varying curvature scale 
(7?o), whereas the distance to the last scattering Xvtc is constant. 
Hence we plot the likelihood as a function of ;^rec/^o. inversely 
proportional to the scale of the fundamental domain. 

Two mulitply-connected realizations of the sky were tested: 
one corresponding to the space in which the last scattering 
sphere can be just inscribed into the fundamental domain, Xvtc = 
Ki, when just the first large angle correlations appear, and the 
second drawn from a somewhat smaller space for which Xmc - 
Rq. We see detections in both cases, stronger as the fundamental 
domain shrinks relative to;^rec- We also calculate the likelihood 
for a model known to be simply-connected. Note that the likeli- 
hood in compact models generically shows a slight increase rela- 
tive to that for the limiting simply connected space as one brings 
the size of the fundamental domain down to the size of the last 
scattering surface (Xrcc ~ "^i), followed, in the absence of sig- 
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Fig. 11: Test for likelihood detectability of compactified space 
for the example of a dodecahedral (/*) closed universe. The 
vertical axis shows the log-likelihood relative to the largest 
model considered. Diff'erent size models are tested against two 
HEALPix A^side = 16 temperature realizations drawn from the 
model with;^rec/^o = 0.314 = (blue) and;^rec/^o = 0.361 
(black). No noise is added and the common mask has been ap- 
plied. Both cases show detection relative to the likelihoods com- 
puted for the isotropic sky realization drawn from the fiducial flat 
infinite universe (red) and compact models of wrong sizes. The 
detection is stronger for the smaller space. Dots mark the posi- 
tions of the models for which the likelihoods were computed. 
Values are given for the orientations of the models which maxi- 
mize the likelihood. The vertical lines show characteristic scales 
of the fundamental domain of the models in the units of curva- 
ture, from smaller to larger, "Ri/Ro, Rm/Ro and Ru/Ro- The vari- 
able XrQc/Ro gives the size of the last scattering surface in the 
same units. The Rq ^ oo limit corresponds to the flat simply- 
connected space. 



nal in the map, by a rapid drop as soon as the models smaller 
than Xyqc are applied. This small increase is also present in the 
fiducial exactly isotropic sky, a single realization of which is 
shown in the figure, but is a generic feature irrespective of the 
topology being tested (occurring also in models with R < Ri), 
and thus should not be taken as an indication for compact topol- 
ogy. The reason for the increase is the possibility of aligning the 
model with a weak anisotropic correlation feature with chance 
patterns of a single sky realization. However the fit drastically 
worsens as soon as the correlation features in a model become 
pronounced. Moreover, the feature becomes considerably less 
significant when the likelihood is marginalized over the orienta- 
tion (Euler angles) of the fundamental domain. 

All of these results (KL divergences and likelihoods) were 
computed with ^^ax = 40, corresponding approximately to 
A^side = 16, indicating that this is more than adequate for de- 
tecting even relatively small fundamental domains such as the 
r[2, 2, 2] case simulated above. We also calculate Jkl between 
the correlation matrices for the T [1,1,1] torus (as a proxy for 
the simply-connected case) and the r[5,5,5] torus, as a func- 
tion of the maximum multipole ^max used in the calculation 
of the correlation matrix: we find that Jkl continues to in- 
crease beyond ^max = 60. Thus, higher-resolution maps (as 
used by the matched-circles methods) contain more informa- 
tion, but with the very low level of noise in the Planck CMB 
maps, ^max = 40 would nonetheless give a robust detection of 
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a multiply-connected topology, even with the conservative fore- 
ground masking we apply. 

We note that it is difficult to compress the content of these 
likelihood figures down to limits upon the size of the funda- 
mental domain. This arises because it is difficult to provide a 
physically-motivated prior distribution for quantities related to 
the size of the fundamental domain. Most naive priors would 
diverge toward arbitrarily large fundamental domain sizes or 
would otherwise depend on arbitrary limits to the topological 
parameters. 

5.3.2. Bianchi 

The ANICOSMO code (McEwen et al. 2013) is used to perform a 
Bayesian analysis of Bianchi VII/^ models, which has been ex- 
tensively validated by McEwen et al. (2013) already; we briefly 
summarise the validation performed for the masked analysis. In 
McEwen et al. (2013) a CMB map is simulated, in which a sim- 
ulated Bianchi temperature map with a large vorticity (i.e., am- 
plitude) is embedded, before applying a beam, adding isotropic 
noise and applying a mask. Both the underlying cosmological 
and Bianchi parameters used to generate the simulations are 
well recovered. For this simulation the coupled Bianchi model 
is favoured over ACDM, with a log-Bayes factor of A In £^ ~ 50. 
As expected, one finds that the log-Bayes factor favours ACDM 
in simulations where no Bianchi component is added. For further 
details see McEwen et al. (2013). 



6. Results 

We now discuss the results of applying the circles-in-the-sky and 
likelihood methods to Planck data to study topology and Bianchi 
VII/j cosmologies. 

6.1. Topology 

Neither the likelihood method nor the circles-in-the-sky search 
find evidence for a multiply-connected topology. We present the 
likelihood for various models. In Fig. 12 we show the likeli- 
hood (marginalized over amplitude and maximized over orien- 
tation of the fundamental domain) for the cubic torus, fixing the 
background cosmology to the best-fit flat Universe Planck model 
(Planck Collaboration XVI 2013). We see that this is maximized 
for L > 2x rec? i-G-? showing no evidence for non-trivial topol- 
ogy. Note that the likelihood shows mild features as the size 
goes through the other scales associated with the topology, in 
particular a small increase in the likelihood when the scale of 
the inscribed sphere is crossed. However, the same increase 
is found when the toroidal model is compared to a single real- 
ization of a strictly isotropic fiducial sky, and thus, should not 
be interpreted as a detection of multi-connected topology. The 
origin of this likelihood behaviour at best fit angles is that the 
freedom of orientation can be used to align small enhancements 
in large-angle correlations in the anisotropic L ^ model with 
random features in the given single realization of the sky. When 
marginalized over all possible orientations the eff'ect is signifi- 
cantly reduced; the slight rise is AlnX ^1.9 from a likelihood 
of P = 650, which is comparable to the numerical noise inherent 
in our stochastic integration. For even smaller spaces, more ex- 
tensive correlations of the temperature can no longer be accom- 
modated and for L < TR^ the likelihood of the T3 cubic toroidal 
model drops quickly, although not strictly monotonically. 



In Figs. 13 and 14 we show the likelihood for the T\L, L, 7] 
chimney and T[L,7,7] slab topologies, which are also maxi- 
mized in the simply-connected limit. The T2 chimney, with only 
two compact dimensions, is less constrained than the T3 cube, 
and the Tl slab, with one compact dimension, even less so. 

We find similar limits for the topologies allowed in a closed 
universe with a locally spherical geometry. In Fig. 15 we show 
the likelihood for the dodecahedral fundamental domain, in 
Fig. 16 for the truncated cube, and in Fig. 17 for the octahe- 
dron. In this case, we do not fix the background cosmologi- 
cal model, but rather account for the geometrical degeneracy 
line which links Hq and Da with Q^. The degeneracy rela- 
tions are approximated as Qa = 0.691 -h 2.705Q/^ and = 
67.8-h388Il/5:-Fl200D^. As in the toroidal case, there is no detec- 
tion of a small space at the level expected from the simulations 
of Sect. 5. Fundamental domains larger than the last scattering 
diameter are preferred for the dodecahedral and truncated cube 
spaces with somewhat weaker restriction for the octahedral case. 
Note that an observationally motivated prior on Hq or Q.k would 
be yet more restrictive on the fundamental domain size. For all 
three topologies, again as in the toroidal case, the maximum of 
the likelihood at best fit orientation is detected for the finite vol- 
ume spaces with^rec ^ Hi Sit the level Aln£ ^ -\-4 relative to the 
fiducial flat simply-connected model. Since this feature is seen 
in the isotropic fiducial sky as well, we cannot take it as an indi- 
cation of a detection of a multi-connected space. In the case of 
curved spaces we see that this mild increase disappears when we 
consider the likelihood marginalized over orientations. 

We present numerical limits for these flat and positively 
curved spaces in Table 2. Because of the one-sided nature of 
these limits, we characterize the shape of the likelihood by the 
steepness of its fall from the value as the scale of the funda- 
mental domain goes to infinity (i.e., the simply-connected limit). 
Hence, we show limits for AlnX < -5, (roughly equivalent to 
a 3(T — 99 % confidence limit — fall for a Gaussian; because 
of the very steep gradient, the 2cr limits are very similar) and 
AlnX < -12.5 (5(t). Note that the limits diff'er depending on 
whether we marginalize or maximize the likelihood over the ori- 
entation angles. We show lower limits on the quantity "Ri (L/2 
for a torus with edge length L) in units of the last scattering dis- 
tance Xrec (in conventional units, Xrec ~ 14 Gpc for the fiducial 
Planck parameters; Planck Collaboration XVI 2013). In most 
cases, the limits are roughly Ri > Xrcc — the scale of the fun- 
damental domain must be greater than that of the last scattering 
surface. We place the most restrictive limits on the dodecahedron 
with Ri > 1.03;^rec using marginalized values for Aln£ < -5. 
Conversely, the chimney and slab spaces are less constrained as 
the expected correlations are weaker in one or two directions; for 
the slab space, we only constrain Ri = L/2 > 0.5;^rec- 

In Fig. 18 we show the likelihood for the two hyperbolic 
models listed in Table 1, which also show no detection of the 
multi-connected topology. In the hyperbolic case we space the 
range of space sizes by varying Qk while keeping Qa - as 
well as Hq constant at fiducial values. 

All of these results show at least some increase in the like- 
lihood for certain orientations when one of the characteristic 
scales of the fundamental domain (R^, or Ri) just exceed 
the surface of last scattering, and so no longer produces matched 
patterns, but induces extra correlations at large angular separa- 
tions. Chance patterns can then mimic these correlations, and 
this is exacerbated by our conservative sky masks, which allow 
arbitrary patterns in the masked regions. 

The circles-in-the-sky search also show no evidence of a 
multiply connected universe. We show the matched circle statis- 
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Fig. 12: Top: the likelihood as a function of the length of an 
edge of the fundamental domain L for a cubic-torus topology. 
In this figure, Xvtc gives the distance to the surface of recom- 
bination. The data are component-separated CMB temperature 
maps degraded to HEALPix A^side = 16 resolution and smoothed 
with an FWHM = 660' Gaussian filter. The common mask 
of /sky = 0-78 is used. The likelihood is marginalized over 
the amplitude of fluctuations, but maximized over the orienta- 
tion of the fundamental domain. Lines for difl'erent estimates of 
the CMB temperature from Planck data are black: SMICA; ma- 
genta: SEVEM; green: Commander -Ruler; blue: NILC. The red 
line is for a simulated isotropic sky from a fiducial flat simply- 
connected model. Noise has been accounted for but is negligible 
at A^side = 16. The likelihoods are normalized to match the like- 
lihood obtained with the common mask in the Rq ^ oo isotropic 
flat limit. The vertical lines mark the positions where Xrec is 
equal to the characteristic sizes of the fundamental domain, from 
left to right, = L/2, = ^^2L/2 and K = V3L/2. Dots, su- 
perimposed onto the SMICA curve, designate the discrete set of 
models studied. Bottom: zoom into the transitional region near 
XrQc ~ ^i- Black Planck SMICA and red fiducial curves are the 
same as in the top panel. The grey curve (open circles) is the like- 
lihood marginalized over the orientations for the Planck SMICA 
map. Only Hi and are within the scale range shown. 



tic in Fig. 19. We do not find any statistically significant corre- 
lation of circle pairs in any map. As seen in Fig. 6, the mini- 
mum radius at which the peaks expected for the matching statis- 
tic are larger than the false detection level is Ofmin ~ 20°. Thus, 
we can exclude at the confidence level of 99 % any topology 
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Fig. 13: Same as Fig. 12, but for a toroidal space with one large 
dimension fixed at IH^^ and two short dimensions of equal size 
L (approximating the "chimney" space), and "Rj^i are marked 
while = oo 
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Fig. 14: Same as Fig. 12, but for a toroidal space with two large 
dimensions fixed at 7H~^ and one short dimension of variable 
L (approximating the "slab" space), is marked while = 
K = oo. 



that predicts matching pairs of back-to-back circles larger than 
this radius, assuming that relative orientation of the fundamen- 
tal domain and mask allows its detection. This implies that in a 
flat universe described otherwise by the Planck fiducial ACDM 
model, a 99 % confidence-limit lower bound on the size of the 
fundamental domain is L/2 > Xvqc ^^^((^mm) = ^-^^Xrcc = 
13.2 Gpc. This is comparable to the cubic torus limits from the 
marginalized likelihood ratios above (L/2 > 0.92;^rec)- 

6.2. Bianchi 

Masked Planck data are analysed for evidence of a Bianchi Yllh 
component, where the prior parameter ranges adopted are the 
same as those specified by McEwen et al. (2013). The analy- 
sis is performed on the SMICA component- separated map, us- 
ing the mask defined for this method, and is repeated on the 
SEVEM component-separated map for validation purposes (using 
the mask defined for the SEVEM method). The Bayes factors for 
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Fig. 15: Top: the likelihood as a function of the distance to 
last scattering surface in curvature units for a locally spheri- 
cal multiply-connected universe with a dodecahedral (/*) fun- 
damental domain with = 0.317?o. Lines are for different esti- 
mates of the CMB temperature from Planck data as in Fig. 12. 
In this figure, the Xrcc/Ro parameterizes the position of the 
model on the geometrical degeneracy line which links Hq and 
D.A with flK- The degeneracy relations are approximated as 
Qa = 0.691 + 2.705Q^ and Hq = 67.8 -h 388Q^ -h 1200D|. 
The red reference curve is for the random isotropic realization 
from a fiducial flat model. Vertical lines mark when;^i-ec equals 
each of HuUm, and the characteristic scales of the funda- 
mental domain. Bottom: zoom into the transitional region near 
XvQc ~ ^i- Both the likelihood at the best orientation of the do- 
main versus the sky (black for the Planck SMICA CMB map and 
red for the fiducial realization, as in the top panel) and the like- 
lihood marginalized over the orientations for Planck SMICA map 
(gray curve, open circles) are shown. 



the various Bianchi VII/^ models and the equivalent standard cos- 
mological models are shown in Table 3. 

For the phenomenological flat-decoupled-Bianchi model, 
evidence in support of a left-handed Bianchi template is found. 
On the Jeflreys scale (Jeflreys 1961), evidence for this model 
would be referred to as strong for the SMICA map and signifi- 
cant for the SEVEM map. For both SMICA and SEVEM component- 
separated data, recovered posterior distributions for the flat- 
decoupled-Bianchi model are shown in Fig. 20a, where similar 
posterior distributions are recovered for both component sepa- 
ration methods. Recall that the Bianchi parameters are decou- 
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Fig. 16: Likelihood for a constant positive curvature multiply- 
connected universe with a truncated cube (O*) fundamental do- 
main with Ri - 0.397?o. Notation is the same as in Fig. 15. 

Table 2: Lower limits on the size of the fundamental domain for 
difl'erent multiply-connected spaces, in units of the distance to 
the last scattering surface, ;^rec- For the torus, slab, and chimney, 
we present limits on the quantity L/2; in curved spaces, limits 
are on the inscribed- sphere topology scale %i. For the columns 
labelled "max", we maximize the probability over the orientation 
of the fundamental domain; for "marg", we marginalize over ori- 
entation. 



Space 


Quantity 




AlnX < -5 


AlnX< 


-12.5 








max 


marg 


max 


marg 


T3 Cubic Torus 




> 


0.83 


0.92 


0.76 


0.83 


T2 Chimney 


^/(2/\rrec) 


> 


0.71 


0.71 


0.63 


0.67 


Tl Slab 


^/(2;t'rec) 


> 


0.50 


0.50 






Dodecahedron 


^i/;rrec 


> 


1.01 


1.03 


1.00 


1.01 


Truncated Cube 




> 


0.95 


1.00 


0.81 


0.97 


Octahedron 


ndx..c 


> 


0.87 


0.89 


0.87 


0.88 



pled from the standard cosmology in the flat-decoupled-Bianchi 
model, hence for this model and are specific to the 
Bianchi model and should not be compared with standard val- 
ues. The maximum a posteriori (MAP) best-fit template found 
for SMICA component- separated data is shown in Fig. 21b, with 
the diff'erence between this template and the template found in 
WMAP 9-year data (McEwen et al. 2013) shown in Fig. 22. Note 
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Fig. 17: Likelihood for a constant positive curvature multiply- 
connected universe with an octahedral (T*) fundamental domain 
with Ri = OASRq. Notation is the same as in Fig. 15. 

Table 3: Log-Bayes factor relative to equivalent ACDM model 
(positive favours Bianchi model). 



Model 


SMICA 


SEVEM 


Flat-decoupled-Bianchi (left-handed) 


2.8 + 0.1 


1.5 + 0.1 


Flat-decoupled-Bianchi (right-handed) 


0.5 + 0.1 


0.5 + 0.1 


Open-coupled-Bianchi (left-handed) 


0.0 + 0.1 


0.0 + 0.1 


Open-coupled-Bianchi (right-handed) 


-0.4 + 0.1 


-0.4 + 0.1 



that the template found in Planck data is very similar to the tem- 
plate found in WMAP 9-year data (McEwen et al. 2013), which 
in turn is similar to the template first found by Jaff'e et al. (2005). 
However, the template found in WMAP 9-year data (McEwen 
et al. 2013) is only significant in full-sky data, but not when the 
9-year KQ75 WMAP mask (Bennett et al. 2012) is applied. Since 
the Planck SMICA and SEVEM masks are less conservative than 
the KQ75 mask, these findings suggest data near the Galactic 
plane may be playing a considerable role in supporting a Bianchi 
component in Planck data. The SMICA CMB map and a Bianchi- 
sub traded version of this map are also shown in Fig. 21. The 
best-fit parameters of the templates found in Planck SMICA and 
SEVEM component- separated data are displayed in Table 4, for 
both the MAP and mean-posterior estimates. The analysis was 
also performed on a SMICA component-separated Gaussian sim- 




0.50 0.60 0.70 0.80 0.90 1.00 1.10 
Zrec^Rj 




Fig. 18: Likelihood for two constant negative curvature 
multiply-connected universe, top: m004(-5,l); bottom: 
v3543(2,3). Notation is as in Fig. 15 except that only "Ri/Ro is 
shown by vertical lines. 



ulation, yielding a null detection (i.e., no evidence for a Bianchi 
component), as expected. 

For the most physically motivated open-coupled-Bianchi 
model where the Bianchi Yllh model is coupled to the stan- 
dard cosmology, there is no evidence in support of a Bianchi 
contribution. Recovered posterior distributions for the open- 
coupled-Bianchi model are shown in Fig. 20b for both SMICA 
and SEVEM component-separated data. Although the cosmolog- 
ical Bianchi parameters agree reasonably well between these 
diff'erent component-separated data, the posterior distributions 
recovered for the Euler angles diff'er. For SEVEM data, an ad- 
ditional mode of the posterior distribution is found; the mode 
found with SMICA data is still present in SEVEM data but is not 
dominant. Consequently, the best-fit estimates for the Euler an- 
gles diff'er between the SMICA and SEVEM component-separated 
data. Note that the additional mode found in SEVEM data is also 
present in WMAP 9-year data (McEwen et al. 2013). The result- 
ing best-fit parameters for the the open-coupled-Bianchi model 
are displayed in Table 5, while the corresponding MAP best-fit 
maps are shown in Fig. 23. Nevertheless, for both SMICA and 
SEVEM data the Bayes factors computed (Table 3) do not favour 
the inclusion of any Bianchi component for the open-coupled- 
Bianchi model. Planck data thus do not provide evidence in sup- 
port of Bianchi VII/^ cosmologies. However, neither is it possi- 
ble to conclusively discount Bianchi Yllh cosmologies in favour 
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Table 4: Parameters recovered for left-handed flat-decoupled-Bianchi model. Planck data favour the inclusion of a Bianchi compo- 
nent in this phenomenological model. 



Bianchi Parameter 




SMICA 




SEVEM 




A/TAD 

MAr 


Mean 


A/TAD 

MAr 


Mean 




0.38 


0.32 + 0.12 


0.35 


0.31+0.15 




0.20 


0.31+0.20 


0.22 


0.30 + 0.20 


X 


0.63 


0.67 + 0.16 


0.66 


0.62 + 0.23 




8.8 X 10-10 


(7.1 + 1.9) X 10-10 


9.4 X 10-10 


(5.9 + 2.4) X 10-10 


a 


38.8° 


51.3° +47.9° 


40.5° 


77.4° + 80.3° 


p 


28.2° 


33.7° + 19.7° 


28.4° 


45.6° + 32.7° 


y 


309.2° 


292.2° +51.9° 


317.0° 


271.5° + 80.7° 



of ACDM cosmologies. The constraints (oj/H)o < 7.6 x 10"^^ 
(95% confidence level) on the vorticity of the physical cou- 
pled Bianchi VII/, left-handed models and (oj/H)o < 8.1 x 10"^^ 
(95% confidence level) for right-handed models are recovered 
from SMICA component-separated data. It will be informative to 
analyse polarised Planck observations for evidence of Bianchi 
models, which may provide a more definitive conclusion. 



7. Discussion 

We have used the Planck temperature anisotropy maps to probe 
the large-scale structure of spacetime. We have calculated the 
Bayesian likelihood for specific topological models in universes 
with locally flat, hyperbolic and spherical geometries, all of 
which find no evidence for a multiply-connected topology with 
a fundamental domain within the last scattering surface. After 
calibration on simulations, direct searches for matching circles 
resulting from the intersection of the fundamental topological 
domain with the surface of last scattering also give a null result 
at high confidence. These results use conservative masks of the 
sky, unlike previous WMAP results, which used full- sky inter- 
nal linear combination maps (not originally intended for cosmo- 
logical studies) or less conservative foreground masks. Hence, 
the results presented here, while corroborating the previous non- 
detections, use a single, self-consistent, and conservative dataset. 
The masked sky also increases the possibility of chance patterns 
in the actual sky mimicking the correlations expected for topolo- 
gies with a characteristic scale near that of the last scattering 
surface. 

Depending on the shape of the fundamental domain, we find 
^ Xyqc (Table 2) with detailed 99 % confidence limits (con- 
sidering the likelihood marginalized over the orientation of the 
fundamental domain) varying from 0.9xyqc for the cubic torus in 
a flat universe to 1.03xyqc for the dodecahedron in a positively 
curved universe, with somewhat weaker constraints for poorly- 
proportioned spaces that are considerably larger along some di- 
rections. 

Note that these results make use of the expected pixel-space 
correlations as a unique signal of non-trivial topology. Hence, 
although a small fundamental domain will suppress power on 
the largest scales of the CMB, observation of such low power on 
large scales as observed by COBE (Bond et al. 2000), and con- 
firmed by WMAP (Luminet et al. 2003), is not suflftcient for the 
detection of topology. Conversely, because our methods search 
directly for these correlations (and indeed marginalize over the 
amplitude of fluctuations), a slight modification of the back- 
ground FRW cosmology by lowering power in some or all multi- 



poles (Planck Collaboration XV 2013) will not aff'ect the ability 
to detect the correlations induced by such topologies. 

Similarly, using a Bayesian analysis we find no evidence 
for a physical, anisotropic Bianchi VII/^ universe. However, 
Planck data do provide evidence supporting a phenomenological 
Bianchi VII/^ component, where the parameters of the Bianchi 
component are decoupled from standard cosmology. The result- 
ing best-fit Bianchi Yllh template found in Planck data is sim- 
ilar to that found in WMAP data previously (Jaff'e et al. 2005; 
McEwen et al. 2013). However, although this Bianchi compo- 
nent can produce some of the (possibly anisotropic) tempera- 
ture patterns seen on the largest angular scales (see also Planck 
Collaboration XXIII 2013), there is no set of cosmological pa- 
rameters which can simultaneously produce these patterns and 
the observed anisotropics on other scales. Moreover, the param- 
eters of the best-fit Bianchi Yllh template in the decoupled set- 
ting are in strong disagreement with other measurements of the 
cosmological parameters. 

These results are expected from previous measurements 
from COBE and WMAP, but Planck's higher sensitivity and 
lower level of foreground contamination provides further con- 
firmation. We have shown that the results are insensitive to the 
details of the preparation of the temperature maps (in particu- 
lar, the method by which the cosmological signal is separated 
from astrophysical foreground contamination). Future Planck 
measurement of CMB polarization will allow us to further test 
models of anisotropic geometries and non-trivial topologies and 
may provide more definitive conclusions, for example allowing 
us to moderately extend the sensitivity to large-scale topology 
(Bielewicz et al. 2012). 
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tion of circle radius a for the Planck CMB maps estimated us- 
ing Coimnander -Ruler (red short dashed line), NILC (blue long 
dashed line), SEVEM (green dot-dashed line) and SMICA (orange 
three dots-dashed line). Dotted line shows the false detection 
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simulations of the CMB map, smoothed and masked in the same 
way as the data, would yield a false event. The peak at 90° cor- 
responds to a match between two copies of the same circle of 
radius 90° centered around two antipodal points. 
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Fig. 20: Posterior distributions of Bianchi parameters recovered from Planck SMICA (solid curves) and SEVEM (dashed curves) 
component- separated data for left-handed models. Planck data provide evidence in support of a Bianchi component in the phe- 
nomenological flat-decoupled-Bianchi model (panel a) but not in the physical open-coupled-Bianchi model (panel b). 
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Fig. 21: Best-fit template of left-handed flat-decoupled-Bianchi VII/^ model subtracted from Planck SMICA component-separated 
data. Before subtraction, the peak-to-peak variation is ±594 yuK, reduced to ±564 yuK after subtraction. 
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